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NON-PERIODIC SYSTEMS WITH 
CONTINUOUS DIFFRACTION MEASURES 


MICHAEL BAAKE, MATTHIAS BIRKNER, AND UWE GRIMM 

Abstract. The present state of mathematical diffraction theory for systems with con¬ 
tinuous spectral components is reviewed and extended. We begin with a discussion of 
various characteristic examples with singular or absolutely continuous diffraction, and 
then continue with a more general exposition of a systematic approach via stationary sto¬ 
chastic point processes. Here, the intensity measure of the Palm measure takes the role 
of the autocorrelation measure in the traditional approach. We furthermore introduce a 
‘Palm-type’ measure for general complex-valued random measures that are stationary and 
ergodic, and relate its intensity measure to the autocorrelation measure. 


1. Introduction 

The (mathematical or kinematic) diffraction theory of systems in Euclidean space with 
pure point spectrum is rather well understood. Ultimately, this is due to the availability 
of Poisson’s summation formula and its generalisations to the setting of measures (or to 
tempered distributions); see [121 Sec. 9.2] for a systematic exposition. Beyond results on the 
spectral nature, this often also provides explicit formulas for the diffraction measure, such 
as in the cases of lattice-periodic systems and model sets. For these systems, there is also a 
well-understood connection with the Halmos-von Neumann theorem for the corresponding 
pure point dynamical spectrum; see [55l [161 EQl lUl 120] for details as well as [12] and 
references therein for general background. 

As soon as one enters the realm of systems with continuous diffraction spectra (or at least 
with continuous spectral components), the situation changes drastically. As in the case of 
Schrodinger operator spectra [31], much less is known about the plethora of possibilities, 
and there rarely are explicit formulas for the diffraction measures of specific examples. Until 
recently, explicit results were restricted to simple systems of Bernoulli type (hence with 
disorder that leads to independent random variables) or to some paradigmatic examples 
in one dimension (and product systems built from them). 

There has now been some progress towards explicitly computable examples in various 
directions [62l H] [15] . In particular, both for singular and for absolutely continuous cases, 
constructive approaches have been more successful than previously anticipated; compare 
[121 Ch. 10]. Consequently, there is some hope that more systems can be understood 

in this way. This view is also supported by the recent progress in the understanding of 
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the connection between the dynamical and the diffraction spectrum in this more general 
situation; see |20] and references therein. At the same time, such examples will improve 
our intuition about systems with continuous diffraction. Below, this will be reflected by 
several short sketches of characteristic examples (which are covered in more detail in jl2j). 
before we embark on a more systematic setting via general point process theory. Our focus 
is on systems in R'^, which is the primary situation to understand, particularly from the 
applications point of view. Extensions to more general locally compact Abelian groups are 
possible, but will not be discussed here. 


2. Diferaction Measures — a Brief Reminder 

Let a; be a locally hnite (and possibly complex) measure on R*^, which we primarily view 
as a linear functional on the space C'c(R‘^) of continuous functions with compact support 
on R'^, together with some mild extra conditions. In favourable cases, u will be translation 
bounded. By the classic Riesz-Markov representation theorem, we may identify the mea¬ 
sures defined by this approach with regular Borel measures; for a systematic exposition, we 
refer to [HI [16] as well as [121 Chs. 8 and 9] and references therein. Particularly important 
examples comprise the Dirac measure S^:, defined by Sx{g) ■= g{x) for g G Cc(R'^), and 
measures of the form 

( 2 . 1 ) h 

x£S 

which are known as Dirac combs, where S' C R*^ is uniformly discrete. More generally, we 
will also consider objects of the form which can be a measure for a general 

countable set S', then under suitable conditions on the weight function w. Such measures 
are referred to as weighted Dirac combs. 

Recall from [H] or na that, if a; is a measure on R"^, the (inverted-conjugate) measure 
uj is defined by uj{g) := oj(g) for g G C'c(K'^), where 'g{x) := g{—x). Given a measure u, 
consider its autocorrelation measure 


( 2 . 2 ) 


7 = 7,.^:='^®'^, 


where ® denotes the volume averaged (or Eberlein) convolution. The latter is dehned by 


Ul ®UJ 


lim 

r—)-CO 


(jJr * OJr 

Vol(i?r(0)) 


with Br{0) the (open) ball of radius r around the origin and oJr ■= i^Ib^(o)- this stage, 
we assume the existence of the limit. This will be discussed in more detail later. 

If (as in many of our examples) ca is a Dirac comb with lattice support, also 7 will be 
supported on the same lattice (or a subset of it). Concretely, if 


= w 5^ := ^w(n) 5n, 

n^'L 


OJ 
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with a bounded weight function w say, one hnds 7 = with the positive definite function 
7 : Z —)■ C being defined by 


(2.3) 


rj{m) 


:= lim - 

N^oo 2A^ + 1 


= lim 


N^OQ 2iV + 1 


N 

w{n) w{n — m) 

n=-N 

N 

w{n)w{n + m), 

n=-N 


provided that all limits exist. In our exposition below, this existence will follow by suitable 
applications of Birkhoff’s ergodic theorem, applied to the dynamical system of the shift 
action on the orbit closure of the sequence tc or to a similar type of dynamical system; 
compare [16] for a more general setting. One benefit of this approach will emerge via the 
Herglotz-Bochner theorem [5U] . 

The autocorrelation measure 7 is positive definite (or of positive type) by construction, 
which means that 'y{g*g) > 0 for all g G Oc(M'^). It is thus Fourier transformable [2^, and 
the Fourier transform 7 is a positive measure, called the diffraction measure of cu; see [29] for 
the physics behind this notion, and [H] as well as [T^ Ch. 9] for the mathematical theory. 
Within the framework of kinematic diffraction, it describes the outcome of a scattering 
experiment by quantifying how much intensity is scattered into a given volume of d-space, 
and thus is the central object of our interest. By the Lebesgue decomposition theorem, 
there is a unique splitting 

'1 Fpp T Tsc T Tac 

of the diffraction measure into its pure point part 7 pp, its singular continuous part 7 sc and 
its absolutely continuous part 7 ac, with respect to Lebesgue measure A. The pure point part 
comprises the ‘Bragg peaks’ (of which there are at most countably many, so %p is a sum 
over at most countably many Dirac measures with positive weights), while the absolutely 
continuous part corresponds to the diffuse ‘background’ scattering which is given by a 
locally integrable density relative to A. The singular continuous part is whatever remains 
— if present, it is a measure that gives no weight to single points, but is still concentrated 
to an (uncountable) set of zero Lebesgue measure. 

Measures u which lead to a diffraction 7 = 7 pp are called pure point diffractive; examples 
include lattice-periodic measures and measures based on model sets. These have been 
studied in detail in the context of diffraction of crystals and quasicrystals; see [9] for a 
recent review and [121 Chs. 8 and 9] for a systematic exposition. Here, we are concentrating 
on the other two spectral components, which may also carry important information on 
the (partial) order which is present in the underlying structure. Pure point spectra are 
discussed in detail in [221 ESI EZl ESI ESI ESI EZl ES] ; for related spectral problems in the 
context of Schrodinger operators, we refer to [31] . 
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Figure 1. The distribution function F of the classic middle-thirds Cantor 
measure. The construction of the underlying Cantor set is sketched in the 
inset. 

3. Guiding Examples 

As mentioned above, the understanding of systems with continuous diffraction compo¬ 
nents is less developed than that of pure point diffractive ones. Still, a better intuition 
will emerge from a sample of characteristic examples. It is the purpose of this section to 
provide some of them, while we refer to the literature for further ones [I3l|22llll|9l[l5l[l2]. 

3.1. Thue—Morse Sequences. Let us begin with a classic example from the theory of 
substitution systems that leads to a singular continuous diffraction measure with rather 
different features in comparison with the Cantor measure, the latter being illustrated in 
Figure [H Our example has a long history, which can be extracted from [THl Ell El II] • 
We conhne ourselves to a brief summary of the results, and refer to min] and references 
therein for proofs and details. 

The classic Thue-Morse (TM) sequence can be dehned via the one-sided hxed point 
V = no'^i'^2 • • • (with Uq = 1) of the primitive substitution rule 


Q\ 


1 ^ 11 
1 ^ 11 












CONTINUOUS DIFFRACTION 


5 


on the binary alphabet {1,1}. The hxed point is the limit (in the obvious product topology) 
of the (suitably embedded) iteration sequence 

1 1 -^ ll 1 -^ 1111 1 -^ 11111111 1 -^ ... —y V = g{v) = VqViV 2 V^ ... 

and has a number of distinctive properties [DEH], for instance 

^ ^_j^^sum of the binary digits of i 

• ^ 2 i = Vi and V 2 i+i = for all i G Nq; 

• V = VqV 2 V^ ... and v = ... 

• T is (strongly) cube-free (and hence non-periodic). 

Here, we dehne 1 = 1 and identify 1 with —1, also for the later calculations with Dirac 
combs. A two-sided sequence w can be defined by 

\vi, for i > 0, 
wn) = < 

for i < 0, 

which is a fixed point of because the seed = 1|1 is a legal word (it occurs in p^(l)) 

and w = Q^{w). The (discrete) hull X = Xrpj^ of the TM substitution is the closure of the 
orbit of w under the shift action, which is a subset of {±1}^ and hence a compact space. 
The orbit of any of its members is dense in X. We thus have a topological dynamical 
system (X, Z) that is minimal. When equipped with the standard Borel cr-algebra, the 
system admits a unique shift-invariant probability measure z/, so that the corresponding 
measure theoretic dynamical system (X, Z, v) is strictly ergodic mm- 
Any given w G X is mapped to a signed Dirac comb u via 

oj = w{n) 6n- 

The image of X is a space of translation bounded measures that is compact in the vague 
topology. We inherit strict ergodicity via conjugacy, and thus obtain an autocorrelation of 
the form of Eq. fl2.2p with coefficients T]{m) as in Eq. fl2.3p . In fact, this autocorrelation 
does not depend on the choice of the element from X, so that we may choose the fixed 
point w from above for the concrete analysis. Due to the nature of w, the coefficients can 
alternatively be expressed as 

1 N-l 

r]{m) = lim — y^VnVn+m 

N^OO 7V ' ^ 
n=0 

for m > 0, together with ri{—m) = r]{m). It is clear that //(O) = 1, and the scaling relations 
of V lead to the recursions 

7](2m) = T]{m) and 
ri{2m + l) = —\{ri{m) + ri{m+l)') ^ 


(3.1) 
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which are valid for all m G Z. In particular, the second relation, used with m = 0, implies 
77 ( 1 ) = — which can also be calculated directly. 

Since r] : Z —> C is a positive dehnite function with 77 ( 0 ) = 1, there is a unique 
probability measure /i on the unit circle (which we identify with the unit interval here) 
such that 

(3.2) ,,(m) = [ dM , 

Jo 

which is a consequence of the Herglotz-Bochner theorem [501 Thm. 1.7.6]. Since u is 
supported on Z, the corresponding diffraction measure 7 is 1-periodic, which follows from 
[3l Thm. 1]; see also [121 Sec. 10.3.2]. One then finds the relation 

7 = /i * 


with the measure /i from Eq. fl3.2l) . appropriately interpreted as a measure on [0,1) and 
hence also on M. Clearly, one also has /^ = 7 |[o i)- One can now analyse the spectral type 
of 7 via that of the finite measure /i, where we follow m-, see also [ 68 l | 20 ] . 

Dehning S{N) = X]m=-Ar(^(’^)) ’ ^ two-step calculation with the recursion (ED estab¬ 
lishes the inequality E{4:N) < ^E{2N) for all iV G N. This implies limN^oo ^{N)/N = 0, 
wherefore Wiener’s criterion [79l|50], see also [121 Prop. 8.9], tells us that /i is a continuous 
measure, so that 7 cannot have any pure point component. Note that the absence of the 
‘trivial’ pure point component of 7 on Z is due to the use of balanced weights, in the sense 
that 1 and —1 are equally frequent. Consequently, the average weight is zero, and the 
claim follows from [121 Prop. 9.2]. 

Let us now dehne the distribution function F by F{x) = /i([0,x]) for x G [0, 1], which 
is a continuous function that defines a Riemann-Stieltjes measure [53l Ch. X], so that 
dF = fi. The recursion relation for 77 now implies na the two functional relations 

dF(|)±dF(7i) = dF(a:), 


which have to be satished by the ac and sc parts of F separately, because we have jj,^^ T 
in the measure-theoretic sense; see [72l Thm. 1.20] or [531 Thm. VII.2.4]. Therefore, defining 

7 ac(m) = l\^~dF,,{x), 


Jo 

we know that the coefficients must satisfy the same recursions fl3.ip as r]{m), 

possibly with a different initial condition 773 ^( 0 ). The classic Riemann-Lebesgue lemma 
[5U1 Thm.1.2. 8 ] states that limm^.±oo hac(’R) = 0. But this limit is only compatible with 
hac(O) = 0, because 77ac(l) = — |7ac(0) and rjad^m) = rjadm) for all 777 G M, so that we 
must have r]ac = 0. This means Fac = 0 by the Fourier uniqueness theorem, wherefore p 
and hence 7 (neither of which is the zero measure) are purely singular continuous. The 
resulting distribution function F is illustrated in Figure [21 Note that F can consistently 
be extended to a continuous function on M via F{x + n) = F{x) + n for n G Z and then 
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Figure 2. The strictly increasing distribution function of the classic, purely 
singular continuous TM measure on [0,1]. 

defines 7 via dF = 7 in the Lebesgue-Stieltjes sense. The function F can efficiently be 
calculated by means of the uniformly converging Volterra iteration 

1 

(3-3) Fn+ii.x) = 2 Jo “ cos( 7 ri/))F„'(i/) dy 

with Fo{x) = X. In contrast to the Devil’s staircase of Figure [H the TM distribution 
function is strictly increasing, which means that there is no plateau (which would indicate 
a gap in the support of 7 ); see [3112] and references therein for details and further properties 
of F. So far, we have obtained the following result. 

Theorem 3.1. Let w be any element of the Thue-Morse hull X = X-pj^, the latter rep¬ 
resented as a closed subshift of {±1}^, and consider the eorresponding Dirac comb wS^. 
Then, its autocorrelation 7 exists and is given fey 7 = 7(5^ with rj being defined by Eq. fld.ip 
together with the initial condition 7 ( 0 ) = 1 . 
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The diffraction measure is ^ * 6-^, where fi is the purely singular continuous proba¬ 

bility measure from Eq. fl3.2l) . In particular, 7 is purely singular continuous as well. □ 

To go one step further, Eq. fl3.3p defines an iteration sequence of distribution functions 
for absolutely continuous measures that converges towards the TM measure in the vague 
topology. Writing dFn{x) = fn{x)dx, one hnds 


n—1 

fn{x) = JJ (1 - COs(2™+Vx)), 

m=0 

which, in the vague limit as u —)■ 00 , gives the well-known Riesz product representation 
of the TM measure; compare [M] for details and [80] for general background on Riesz 
products. 

The TM sequence is closely related to the limit-periodic period doubling (pd) sequence, 
compare mm and references therein, via the (continuous) sliding block map dehned by 

(3.4) 0: 11,11^6, 

which results in an exact 2-to-l surjection from the hull to Xp^. The latter is the hull 
of the period doubling substitution dehned by 

(3.5) at—)-a 6 , 6 ha aa. 

Viewed as topological dynamical systems, this means that (Xp^, Z) is a factor of (X-pj^, Z). 
Since both are strictly ergodic, this extends to the corresponding measure-theoretic dy¬ 
namical systems. The period doubling sequence can be described as a regular model set 
with a 2-adic internal space |2H 123] and is thus pure point diffractive. This pairing also 
explains a phenomenon observed in [78], namely that the dynamical spectrum of the TM 
system is richer than its diffraction spectrum. By the dynamical (or von Neumann) spec¬ 
trum, we mean the spectrum of the unitary operator induced by the shift on the Hilbert 
space T^(X, v), where v is the unique shift-invariant probability measure on X; see [HB] for 
more. Here, the pure point part of the dynamical spectrum is the ring Z[|], which is not 
even hnitely generated (and only the ‘trivial’ part Z is detected by the diffraction measure 
of the TM system with general weights). In fact, our above measure /i from Theorem 13.11 
represents the maximal spectral measure in the ortho-complement of the pure point sec¬ 
tor [ 6 HIEQ]. The missing pure point part, however, is fully recovered via the diffraction 
measure of Xp^; see [ 12 ] for details and [ 20 ] for a general discussion of this phenomenon. 

Various generalisations of this result are known by now. First of all, and perhaps not 
surprisiningly, this generalises to an entire family of bijective, binary substitutions [5]. 
Moreover, extensions to higher dimensions are also possible, including the explicit nature 
of the resulting diffraction measure; compare [iQlin] and references therein. 
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3.2. Rudin—Shapiro Sequence. The (binary) Rudin-Shapiro (RS) chain is a bi-in£nite 
deterministic sequence, with polynomial (in fact linear) complexity function and thus zero 
entropy. It can be described recursively as w = (w(n))„gg with w{n) G {±1}, with initial 
conditions tc(—1) = —1, tc(0) = 1 and the recursion 


(3.6) 


ic(4n + i) 


w{n), for i G { 0 , 1 }, 

(— 1 )”+^ w(? 7 ,), for i G { 2 , 3 }, 


which determines w{n) for all n G Z. The orbit closure of w under the shift action is the 
(discrete) RS hull Alternatively, one can start from a primitive substitution on a 

4-letter alphabet (via a i—)■ ac, 6 i—)■ dc, c i—)■ ab and d i—)■ db) and define a quaternary hull, 
which then maps to the binary hull via a simple reduction to two letters (for instance via 
a, c I—)■ 1 and 6, d i—)■ —1); compare [HEH] or [121 Sec. 4.7.1] for details. The two hulls dehne 
topologically conjugate dynamical systems, with local derivation rules in both directions; 
see [121 Rem. 4.11]. 

The shift action on Xj^g is strictly ergodic, so that one can dehne functions r],'d: Z — > C 
via 


rjim) 

'd{m) 


lim - 

N^oo 2N +1 


lim - 

N^OO 2N +1 


N 

w{n) w{n — m) and 

n=-N 

N 

(—l)"w(n) w(n — m), 

n=-N 


where all limits exist due to unique ergodicity (which is best formulated on the level of 
the 4-letter alphabet mentioned above). In particular, one hnds 77 ( 0 ) = 1 and 79(0) = 0. 
The recursive structure of Eq. fl3.6p now implies the validity of a closed set of recursive 
equations mu, namely 


together with 


77(4777) 

_ l+(-l)"* 
2 

■ 77(777), 



77(4777-1-1) 

_ l-(-l)"* 
“ 4 

■ 77(777) -f — 

79(777) — 

179(777 -M), 

77(4777-1-2) 

= 0, 




77(4777-1-3) 

_ l+(-l)"‘ 
4 

+ 

1 

- 79(777] 

) + i79 ( 777 +l). 

79(4777) : 

= 0, 





79(4m-|-l) = — + + 


79(4m-|-2) = 


— ^ 'd{m) + ^'d{m+l), 


79(4m-|-3) = 


i+(-i)^ 


77(777-1-1) — ^ 79(777) -I- j'd{m+l). 


4 
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which hold for all m G Z; see |T 2 l Sec. 10 . 2 ] for details. A careful inspection shows that 
the unique solution of this set of equations, with the initial conditions mentioned above, 
is ri{m) = 5^0 together with 'd(m) = 0 for all m E 1^. Hence, despite the deterministic 
nature of the RS sequence, the autocorrelation measure is simply given by 7 j^g = 5g, so that 
= A, where A again denotes Lebesgue measure. Alternatively, the result also follows 
from the exposition in |6H1 167] . 

Theorem 3.2. Let w be any element of the Rudin-Shapiro hull Xj^g C {±1}^, and con¬ 
sider the corresponding Dirac comb wS^. Then, its autocorrelation exists and is given by 
7rs ~ *^0? diffraction measure = A. □ 

As in the case of the TM sequence, the non-trivial pure point part of the dynamical 
spectrum (which is Z[A] once again) is not ‘seen’ by the diffraction measure, while A 
(with multiplicity 2 ) represents once again the maximal spectral measure in the ortho¬ 
complement of the pure point sector. However, the missing pure point component can be 
recovered by a suitable factor system, the latter obtained via the block map dehned by 
Eq. (13.4p . The corresponding factor is represented by a limit-periodic substitution rule that 
is somewhat reminiscent of the paper folding sequence [1]; see [121 Sec. 10 . 2 ] for a complete 
discussion and [ 20 ] for the general connection between dynamical and diffraction spectra. 
The structure underlying the RS sequence can be generalised to higher-dimensional lattice 
substitutions in a rather systematic way; see |3S] for details. 

3.3. Bernoullisation. Let us begin this discussion by recalling the structure of the full 
Bernoulli shift from the viewpoint of kinematic diffraction. The classic coin tossing process 
leads to the Dirac comb 

u = y^^X{n)6n, 

where the form an i.i.d. family of random variables, each taking values 1 and 

—1 with probabilities p and 1 — p, respectively. By an application of the strong law of 
large numbers (SLLN, see [36] for a favourable formulation), almost every realisation has 
the autocorrelation measure 

7 = (2p - 1)2 + 4p(l - p) 5o , 

and hence (via Fourier transform) the diffraction measure 

7 = (2p - 1)2 dz + 4p(l - p) A . 

Here, we have used the classic Poisson summation formula = 6^; compare [S] and 
references therein, as well as [T 2 [ Sec. 9.2] for a formulation in the diffraction context. 
When p = |, the diffraction boils down to 7 = A. Here, the point part is extinct because 
the average scattering strength vanishes. For proofs, we refer the reader to [ 22111 ], while 
[SI1E2] contain several important and non-trivial generalisations and extensions; see also 
[56] for important related material. 
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The Bernoulli chain has (metric) entropy |28l l35] 

Kp) = -plog(p) - (l-p) log(l-p), 

which is maximal for p = |, with h(|) = log(2). It vanishes for the deterministic limiting 
cases p G {0,1}. For the latter, we have oj = and consequently obtain the pure point 
diffraction measure 7 = 5^, again via Poisson’s summation formula. 

Now, the theory of random variables allows for an interpolation between deterministic 
(binary) sequences and coin tossing sequences as follows. If tc G {±1}^ denotes a deter¬ 
ministic sequence (which we assume to be uniquely ergodic for simplicity), consider the 
random Dirac comb 

(3.7) ujp = '^w{n)X{n)5n, 

where (X(n))^gg is, as above, an i.i.d. family of random variables with values in {± 1 } and 
probabilities p and 1 — p. This ‘Bernoullisation’ of w can be viewed as a ‘model of second 
thoughts’, where the sign of the weight at position n is changed with probability 1 — p; 
compare [121 Sec. 11 . 2 . 2 ]. 

Let w now be the Rudin-Shapiro sequence from above. By a (slightly more complicated) 
application of the SLLN, it can be shown [ 8 ] that the autocorrelation 7 ^ of the Dirac comb 
cjp is then almost surely given by 

Ip = (2p - 1)^ 7rs + 4p(l - P) ^0 = ^0 ) 

irrespective of the value of the parameter p G [0,1]. Recall that two measures with the same 
autocorrelation are called homometric; see [121 Sec. 9.6] for background. Our observation 
thus establishes the following classic result; see [SI El 112] for details. 

Theorem 3.3. The random Dirac combs Up of Equation fl3.7p with real parameter values 
p G [0,1] are {almost surely) homometric, with absolutely continuous diffraction measure 
^ = 7 ^ = A, irrespective of the value of p. In other words, the family [up | p G [0,1]} 
is {almost surely) isospectral. □ 

This result shows that diffraction can be insensitive to entropy, because the family of 
Dirac combs ujp of Eq. fl3.7p continuously interpolates between the deterministic Rudin- 
Shapiro case with zero entropy and the completely random Bernoulli chain with maximal 
entropy log(2). Clearly, the Bernoullisation procedure can be applied to other sequences 
as well, and can be generalised to higher dimensions. For further aspects of entropy versus 
diffraction, we refer to iiiiniiig]. 

3.4. Random Dimers on the Line. Another instructive example [25] is based on certain 
dimer configurations on Z. To formulate it, we follow the exposition in ra and partition Z 
into a close-packed arrangement of ‘dimers’ (pairs of neighbours), without gaps or overlaps. 
Clearly, there are just two possibilities to do so, because the position of the first dimer fixes 
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that of all others. Next, decorate each dimer randomly with either (1,-1) or (—1,1), with 
equal probability. This results in patches such as 

• • • [+ -][- +][- +][+ -][- +][- +][- +][+ -][+ -] • • • 

• • • [- +][+ -][+ -][- +][+ -][+ -][+ -][- +][- +][+ -] • • • 
where the dimer boxes are indicated by brackets. The set of all decorated sequences dehned 
in this way is given by 

X = [w e {±1}^ I M{w) C 2Z or M{w) C 2Z + l} , 

where M{w) := {n G Z | w{n) = w{n + 1)}. Note that M{w) is empty precisely for the 
two periodic sequences that are dehned by w{n) = ±(—1)"' for n G Z. Clearly, X C {±1}^ 
is closed and hence compact. 

Let tc G X and consider the corresponding signed Dirac comb on Z with binary weights 
w{n) G {±1}. One can then show (again via the SLLN) that the corresponding autocor¬ 
relation almost surely exists and is given by [25] 

(3.8) 7 = - 1(5^ + 5_^). 

The corresponding diffraction measure is then 

(3.9) 7 = (l — cos(27rfc)) A , 

which is again purely absolutely continuous. Here, the (smooth) Radon-Nikodym density 
relative to A is written as a function of k. Note that the diffraction measure for general 
weights /i_|_ and h- is given by 

_ \h+ + h_|2 , |h+ - h_|2 _ 

= - 1 -' 5 ^ +- 1 - 7 

with 7 as in Eq. fl3.9|) . In particular, the measure ^ shows only the ‘trivial’ pure point 
diffraction contribution that arises as the consequence of Z being the support of the 
weighted measure under consideration. The same phenomenon also occurs for general 
(non-balanced) TM and RS sequences; compare [121 Rems. 10.3 and 10.5]. 

On first sight, the system looks disordered, with entropy |log(2). This seems (qualita¬ 
tively) reflected by the diffraction. However, the system also dehnes a measure-theoretic 
dynamical system under the action of Z, as generated by the shift. As such, it has a dy¬ 
namical spectrum that does contain a pure point part, with eigenvalues 0 and A we refer 
to [HE] for general background on this concept, and to ESI for the actual calculation of the 
eigenfunctions. The extension to a (continuous) dynamical system Xc under the general 
translation action of R is done via suspension; see [28l Ch. 11.1] (where the suspension is 
called a special flow) or [35] for general background. 

This Ending suggests that some degree of order must be present that is neither visible 
from the entropy calculation nor from the difiraction measure alone. Indeed, in analogy 
with the situation of the TM and the RS sequence, one can define a factor of the system 
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by a sliding block map 0: X —)■ {±1}^ defined by {(j)w){n) = —w{n)w{n + 1). It maps X 
globally 2:1 onto 

Y = 4)(X) = {n G {±1}^ I v{n) = 1 for all n G 2Z or for all n G 2Z + l}. 

The suspension Yc (for the action of M) is defined as above. The mapping 0 extends 
accordingly. 

The autocorrelation and diffraction measures of the signed Dirac comb v6^ for an element 
n G Y are almost surely given by 

7 = 2*^0 + 2*^22 and 7 = -A + - 5^/2 • 

The diffraction of the factor system Y uncovers the ‘hidden’ pure point part of the dy¬ 
namical spectrum, which was absent in the purely absolutely continuous diffraction of the 
signed Dirac comb with re G X. In summary, we have the following situation [2^ [20] . 

Theorem 3.4. The diffraction measure of the close-packed dimer system X with balanced 
weights is purely absolutely continuous and given by Eq. fl3.9p . which holds almost surely 
relative to the natural invariant measure of the system. 

The dynamical spectrum of the continuous close-packed dimer system Xc under the trans¬ 
lation action of M contains the pure point part 'Ll2 together with a countable Lebesgue 
spectrum. 

The non-trivial part L-\-k of the dynamical point spectrum is not reflected by the diffrac¬ 
tion spectrum of X^ but can be recovered via the diffraction spectrum of a suitable factor, 
such as Yc. □ 

As in the case of the Thue-Morse system, where the missing pure point part of the 
dynamical spectrum is recovered by the diffraction measure of the period doubling factor, 
we thus see that and how we can recover the missing eigenvalue via a generalised 2 -point 
function. This observation can be extended to symbolic systems over finite alphabets and 
also to uniquely ergodic Delone dynamical systems of finite local complexity; see [20] for 
details. 

3.5. Ledrappier’s Shift Space. For a long time, people had expected that higher di¬ 
mensions are perhaps more difficult, but not substantially different. This turned out to be 
a false premise though, as can be seen from the now classic monograph [75] . 

In our present context, we pick one characteristic example, the system due to Ledrappier 
[54] , to demonstrate a new phenomenon. We follow the brief exposition in [10] and consider 
a specific subset of the full shift space {± 1 }^^, defined by 

(3.10) Xl = {tc G {±1}^^ I w{x) w{x -\- e^) w{x -|- 63 ) = 1 for all x G 

where and 63 denote the standard Euclidean basis vectors in the plane. On top of being 
a closed subshift, Xl is also an Abelian group (here written multiplicatively), which then 
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comes with a unique, normalised Haar measure. The latter is also shift invariant, and the 
most natural measure to be considered in our context; see also the reformulation in terms 
of Gibbs (or equilibrium) measures in [71]. 

The system is interesting because the number of patches of a given radius (up to trans¬ 
lations) grows exponentially in the radius rather than in the area of the patch. This 
phenomenon is called entropy of rank 1, and indicates a new class of systems in higher 
dimensions. More precisely, along any lattice direction of 7?, the linear subsystems essen¬ 
tially behave like one-dimensional Bernoulli chains. It is thus not too surprising that the 
diffraction measure satishes the following theorem, though its proof [2S] has to take care 
of the special directions connected with the dehning relations of Xl. 

Theorem 3.5. If w is an element of the Ledrappier subshift Xl of Eq. (13.101) . the corre¬ 
sponding weighted Dirac comb has diffraction measure X, which holds almost surely 
relative to the Haar measure of Xl. □ 

So, the Ledrappier system is homometric to the (full) Bernoulli shift on {±1}^^, which 
means that an element of either system almost surely has diffraction measure A. As men¬ 
tioned before, via a suitable product of two Rudin-Shapiro chains, also a deterministic 
system with diffraction A exists. This clearly demonstrates the insensitivity of pair corre¬ 
lations to the (entropic) type of order or disorder in the underlying system; see also |8]. 
Due to the dehning relation in Eq. (13.101) . it is clear that certain three-point correlations in 
the Ledrappier system cannot vanish, and thus make it distinguishable from the Bernoulli 
shift. 

Although correlation functions of third order can resolve the situation in this case (and 
in many other examples as well |32l|57|), one can consider other dynamical systems (such 
as the (x2, x3)-shift [26]) that share almost all correlation functions with the Bernoulli 
shift on [0,1]^^. This is a clear indication that our present understanding of ‘order’ is 
incomplete, and that we still lack a good set of tools for the detection and classihcation of 
order. For a recent alternative based on direct space statistics, we refer to [6]. 

3.6. Random Matrix Ensembles. Another interesting class of random point sets derives 
from the (scaled) eigenvalue distribution of certain random matrix ensembles; see [H] and 
references therein. The global eigenvalue distribution of random orthogonal, unitary or 
symplectic matrix ensembles is known to asymptotically follow the classic semi-circle law. 
More precisely, this law describes the eigenvalue distribution of the underlying ensembles 
of symmetric, Hermitian or (symplectically) self-dual matrices with Gaussian distributed 
entries. The corresponding random matrix ensembles are called GOE, GUE and GSE, 
with attached /^-parameters 1, 2 and 4, respectively. They permit an interpretation as 
a Goulomb gas, where /3 is the power in the central potential; see [21 163] for general 
background and [SH [38] for the results that are relevant here. 
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Figure 3. Absolutely continuous part of the autocorrelation (left) and the 
diffraction (right) for the three random matrix derived point set ensembles on 
the line, with f3 G {1,2,4}. On the left, the oscillatory behaviour increases 
with 13. On the right, (3 = 2 corresponds to the piecewise linear function 
with bends at 0 and ±1, while f3 = A shows a locally integrable singularity 
at ±1. The latter reflects the slowly decaying oscillations on the left. 

For matrices of dimension N, the semi-circle has radius ^J2N/^T and area N. Note that, 
in comparison with [63] , we have rescaled the density by a factor 1 / so that we really 
have a semi-circle (and not a semi-ellipse). To study the local eigenvalue distribution for 
diffraction, we rescale the central region (between ±1, say) by ^J2Nj-K. This leads, in the 
limit as — )■ cxo, to an ensemble of point sets on the line that can be interpreted as a 
stationary, ergodic point process of intensity 1; for (3 = 2, see |2l Ch. 4.2] and references 
therein for details. Since the underlying process is simple (meaning that, almost surely, no 
point is occupied twice), almost all realisations are point sets of density 1 . 

It is possible to calculate the autocorrelation of these processes, on the basis of Dyson’s 
correlation functions [M]- Though these functions originally apply to the circular ensem¬ 
bles, they have been adapted to the other ensembles by Mehta [63] . For all three ensembles 
mentioned above, this leads to an autocorrelation of the form 

(3.11) 7 = ^0 + (1 -/(kl))A 

where / is a locally integrable function that depends on /3; see HI for the explicit formulas, 
and the left panel of Figure [3] for an illustration. 

The diffraction measure is the Fourier transform of 7 , which has also been calculated in 
[341 [63] . Recalling ^0 = -^ and A = fhe result is always of the form 

(3.12) 7 = ^0 + (1 - h{k))\ = ^0 + h{k) A, 

where b = f. The Radon-Nikodym density h depends on (3 and is summarised in [14] . 
Figure [3| illustrates the result for the three ensembles. 
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A similar approach is possible on the basis of the eigenvalues of general complex random 
matrices. This leads to the ensemble studied by Ginibre [63], which is also discussed in 
[H]. One common feature of the resulting point sets is the effectively repulsive behaviour 
of the points, which leads to the ‘dip’ around 0 for 7. For the two systems mentioned in this 
section, we omit the formulation of the full results and refer the reader to for details. 
Further developments around determinantal and related point processes are described in 
reference |15j . 


4. The Renewal Process 

A large and interesting class of processes in one dimension can be described as a renewal 
process jaTiiiiiii]. Here, one starts from a probability measure fi on M+ (the positive real 
line) and considers a machine that moves at constant speed along the real line and drops 
a point on the line with a waiting time that is distributed according to p. Whenever this 
happens, the internal clock is reset and the process resumes. Let us (for simplicity) assume 
that both the velocity of the machine and the expectation value of p are 1 , so that we end 
up with realisations that are, almost surely, point sets in R of density 1 (after we let the 
starting point of the machine move to — 00 , say). 

Clearly, the resulting process is stationary and can thus be analysed by considering all 
realisations which contain the origin. Moreover, there is a clear (distributional) symmetry 
around the origin, so that we can determine the corresponding autocorrelation 7 of almost 
all realisations from studying what happens to the right of 0. Indeed, if we want to 
know the frequency per unit length of the occurrence of two points at distance x (or the 
corresponding density), we need to sum the contributions that x is the hrst point after 
0, the second point, the third, and so on. In other words, we almost surely obtain the 
autocorrelation 

(4.1) 7 = 5o + u + u 

with z/ = /i + /i=t:/i + /i*/i=t:/i + ..., where the proper convergence of the sum of iterated 
convolutions follows from [H Lemma 4] or from [T21 Sec. 11.3]. Note that the point measure 
at 0 simply reflects the fact that the almost sure density of the resulting point set is 1 . 
Indeed, i/ is a translation bounded positive measure, and satishes the renewal relations 
(compare [ST] Ch. XL9] or [H Prop. 1] for a proof) 

(4.2) z/ = fi + fi* v and (1 — p) P = p, 

where p is a uniformly continuous and bounded function on R. The second equation 
emerges from the hrst by Fourier transform, but has been rearranged to highlight the 
relevance of the set S' = {fc | /i(fc) = 1} of singularities. In this setting, the measure 7 of 
Eq. fl4.ll) is both positive and positive dehnite. 

Based on the structure of the support of the underlying probability measure /r, one 
can determine the diffraction of the renewal process explicitly. To do so for a probability 
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measure /r on M+ with mean 1 , we assume the existence of a moment of /i of order 1 + £ 
for some £ > 0; we refer to [1] for details on this condition. The diffraction measure of the 
point set realisations of the stationary renewal process based on /r almost surely is of the 
form 

7 = %p + (1 - h) A, 

where h is a locally integrable function on M that is continuous almost everywhere. The 
pure point part is trivial, meaning 7 = 5q, unless the support of fi is contained in a lattice. 
The details are stated below in Theorem 14.11 Proofs of these claims as well as further 
results can be found in [U [TH [ 12 ] . 

The renewal process is a versatile method to produce interesting point sets on the line. 
These include random tilings with finitely many intervals (which are Delone sets) as well 
as the homogeneous Poisson process on the line (where /i is the exponential distribution 
with mean 1); see [U Sec. 3] for explicit examples and applications. In particular, if one 
employs a suitably normalised version of the Gamma distribution, one can formulate a one- 
parameter family of renewal processes that continuously interpolates between the Poisson 
process (total positional randomness) and the lattice Z (perfect periodic order); compare 
[H Ex. 3] for more. The general result reads as follows. 


Theorem 4.1. Let g he a probability measure on M_|_ with mean 1, and assume that a 
moment of g of order 1 -f e exists for some £ > 0. Then, the point sets obtained from the 
stationary renewal process based on g almost surely have a diffraction measure of the form 


7 = lpp + 0--h)\, 


where h is a locally integrable function on R that is continuous except for at most countably 
many points {namely those of the set S' = {fc | 'g{k) = 1}). On R \ S', the function h is 
given by 


h{k) 


2 (|2(fc)P - Re{g{k))) 
|1 - g{k)\‘^ 


Moreover, the pure point part is 


7pp 



if supp(p) is not a subset of a lattice, 
otherwise, 

where bT> is the coarsest lattice that contains supp(p). 


□ 


In one dimension, the renewal process allows an efficient derivation of the diffraction of 
random tilings, which we briefly summarise now. 


5. Random Tilings 

The deterministic Fibonacci chain can be defined by the primitive substitution rule 
a HA a 6 , 6 HA a, which defines a strictly ergodic (discrete) hull. When a and b are replaced 
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by intervals of length r = |(1 + \/5 ) and 1, respectively, the left endpoints of the intervals 
define a model set (or cut and project set). The corresponding Dirac comb leads to the 
pure point diffraction measure 

^=5^ I{k)6, 

with intensities I{k) = Here, ^ ^ is the density of the point set, and 

k' denotes the algebraic conjugate of fc, which is dehned on the held Q(\/5) by \/5 i—)■ —\/5 
and acts as the A-map for the underlying model set description. In particular, the diffraction 
is the same for all Dirac combs of the Fibonacci hull; see [121 Sec. 9.4.1] and references 
therein for details. An illustration is shown in the upper panel of Figure HI 

The corresponding random tiling ensemble consists of all tilings of the real line by the 
two types of intervals. For a direct comparison, it makes more sense to only consider those 
tilings with the same relative frequency of interval types, which means frequencies 1/r and 
1 /r^ for the long and the short interval, respectively. 

The diffraction of a typical Dirac comb out of this class was originally derived in [13], 
but can also be obtained via an application of the renewal structure from Theorem 14.11 
This leads to 

tS = ^i„ + fcA 

with the Radon-Nikodym density 

h(k'\ = ^ + ^ _ (sin( 7 rfc/r))^ _ 

5 r^(sin(7r/cr))^ + r(sin(7rA:))^ — (sin(7rfc/r))^ 

Except for the trivial Bragg peak at /c = 0, the diffraction measure is thus absolutely 
continuous. Still, the resemblance between this function and the diffraction of the perfect 
Fibonacci chain is remarkable, as can be seen from Figure HI 

The situation in dimensions d > 2 is less favourable from a mathematical perspective, 
although one has a rather clear intuition of what one should expect dal [70], based on 
solid scaling arguments. In dimensions d > 3, a mixed spectrum with pure point and 
absolutely continuous components is conjectured, while d = 2 is the critical dimension 
in the sense that random tilings with non-crystallographic symmetries should display a 
singular continuous component; see [121 Sec. 11.6.2] for an example. 

Unfortunately, only few results have been proved so far. Among them are a rigorous 
treatment of planar random tiling ensembles with crystallographic symmetries (such as 
the lozenge tiling and several relatives, see [I31H51H6]), a group-theoretic approach to one 
of the random tiling hypotheses [69l [70] and a treatment of dense Dirac combs with pure 
point diffraction [THEH] that is needed to understand the pure point part of the random 
tiling diffraction in dimensions d > 3. The remaining questions are still open, though there 
is little doubt that the original analysis from [43] is essentially correct. 








CONTINUOUS DIFFRACTION 




Figure 4. The pure point diffraction measure of the perfect Fibonacci chain 
(upper panel) and the absolutely continuous part of the corresponding ran¬ 
dom tiling (lower panel). Bragg peaks (in the upper picture) are shown as 
lines, where the height is the intensity, while the smooth Radon-Nikodym 
density in the lower picture is truncated at a value of 20 to illustrate the 
spikyness. The central peak (of intensity ^^) is omitted in both diagrams. 
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Let US now leave the realm of explicit examples and turn our attention to a more general 
approach of systems with randomness, formulated with methods from the theory of point 
processes; compare @21 EH Eg S] for related aproaches and results. 

6 . Stochastic Point Processes and the Palm Measure 

In this section, we take the viewpoint of a general shift-invariant random measure and 
relate its realisation-wise diffraction to its second moment measure. As such, this section 
is a complex-valued extension of @1 Sec. 5]. 

Let he a locally finite complex-valued measure on (which means that 

/ijj and are both locally finite signed measures). A short calculation reveals that, 
for / G Cc(M'^, C) of the form f = g + ih with real-valued g and h, the measure p can 
consistently be defined via 

Kf) ■= h(/) = 7^(/-) = /^5 r(/-) - ihs(/-), 

where f{x) = /_(a;) with f_{x) = /(—x). In particular, note that 

JI = i/icj and fi{f) = Jl(f) 

hold as expected. The point here is that, after having dealt with the case of real (or signed) 
measures, the extension to complex measures is canonical and consistent. 

To continue, recall the polar representation of a complex measure from @31 Ch. XIIL16]; 
see also m Prop. 8.3]. Given /i, there is a measurable function a^: —)■ [0,27r) such 

that, for / G Gc(M'^, C), one has 

/ f{x)dfi{x) = [ /(x)e‘“^(^M|/x|(x), 

jRd. jRd 

where |/i| is the total variation measure of /i. This means that |/x| is the smallest non¬ 
negative measure such that |/i(A)| < |/i|(A) for any bounded and measurable A, where 
Ihl ^ I/^5 r| + compare [I2l Sec. 8.5.1] and references therein. 

Let Ai denote the C-vector space of all locally finite, complex-valued measures 0 on 
so 0 G AT means |0(A)| < oo for any bounded Borel set A. A sequence (0n)neN C Ai 
converges vaguely to 0 if 0n(/) —> 0(/) as n —)■ oo for all / G Gc(M'^). The space A4 
is closed in the topology of vague convergence of measures (in fact, AT is a Polish space 
with this topology). We let Um denote the a-algebra of Borel sets of AT. The latter can 
be described as the cr-algebra of subsets of AT generated by the requirement that, for all 
bounded Borel sets A C the mapping 0 i—)■ 0(A) is measurable. 

For each t G let Tt denote the translation operator on M'^, as dehned by the mapping 
X t + X. Clearly, T^Tg = and the inverse of Tt is given by = T_t. For 

functions / on the corresponding translation action is defined via Ttf = / o T_t, so 
that (Ttf){x) = f{x — t). Similarly, for 0 G AT, let := 0 o be the image measure 
under the translation, so that {Tx4>){A) = 0(T_3,(A)) = 0(A — x) for any measurable 
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subset A C and f{y) d{T^(j)){y) = f {x + z) d(j){z) = for 

functions. This means that there is a translation action of on Finally, we also have 
a translation action on P(A4), the probability measures on J\4, via (TxQ)(A) = Q(T_xA) 
for any A G and Q G V{Ai). A set A G Um is called invariant (under translations) 
if T-xA = A for all x G 

A (complex-valued) random measure ^ is a random variable (dehned on some probability 
space {0, J-, vr)) with values in Ai, which formally means that ^: O — > M is an 
measurable function. Its distribution is then Q = 7 i o g V{Ai), i.e. the image measure 
of 71 under <P. We will follow the usual practice in probability theory and not make the 
underlying probability space explicit (a canonical choice can in many cases simply be 
0 = M and = Id^^). We will also usually suppress the dependence of ^ on 6^ G 6 > in 
the notation. Integrals over 0 w.r.t. the probability measure tt will be denoted by E, the 
expectation value. 

<d> is called stationary if its distribution Q satishes TxQ = Q for all x G A stationary 
random measure is called ergodic if the shift-invariant cr-algebra is trivial, which means 
that any invariant A has probability 0 or 1 (more generally, one requires Q{A) G {0,1} 
whenever Q{(T_xA) A A) = 0 for all x G compare [301 Def. 10.3.1 and Prop. 10.3.III]). 


In what follows, we generally assume that 

X ^ is a (possibly) complex-valued, stationary and er- 

^ godic random measure on 

which means that there is a decomposition where both and are signed, 

real-valued, stationary, ergodic random measures on M'^. To verify the last statement note 
that, since for any bounded measurable A C 6 i—)■ ^{A) (= ^{6, A)) G C is measurable, 
also ^sr(A) and ^q(A) are measurable as functions of 6. Consider any shift-invariant 
measurable B C Ad real (Ad real denotes the locally hnite signed measures on R'^), then 
{<P I G B} is shift invariant and measurable as well, so P(^k G 5) G {0,1}, and 
analogously for We further assume that ^ is locally square integrable in the sense that 


( 6 . 2 ) 


E 


(|<f*|(.4))'U (|4.3|(.4))- 


< oo 


for all bounded A C R'^, 


where and |^c-| denote the total variation measures of and respectively. 

In analogy with the real-valued case in [H Sec. 5.2], we dehne the second moment 
measure of via 


(6.3) X A') = E[^(A)^(A')] for bounded A, A' G 

hence, for / G Cc(R'^ x R'^, C), 


X"'(/) = E 


f{x,y)dt{x)d0(y) 
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By the shift invariance of the distribution of iP, we have 

/ = / f{x + z,y + z)dy^‘^\x,y) 

JM.'ixRd, 

for all 2 ; G and hence we can factor out this symmetry to obtain the reduced sec- 

( 2 ) 

ond moment measure The latter is a locally hnite complex-valued measure that is 

characterised by 


(6.4) 


^Xll 

for / G Cc(M‘^ X C). By the shift invariance of Lebesgue measure on we equivalently 
have 

(6.5) [ f{x,y)dfi^‘^\x,y) = ^ ^ ^ 


f{x,y)dy^^\x,y) = f {u + v, u) dy^^l{v) dX{u) 



f{u,u-v) dA(M). 

To prove the existence of one can decompose = /r^^-|-i/r^^ into real and imaginary 
parts and then use the well-known real-valued results (compare [30l Lemma 10.4.Ill]) to 
obtain 

ODtam — /igtj j,g(i -I- 1/^3 

Note that is uniquely defined and is a positive definite measure, since 

= [ [ f{u + v)g{u)dylll{v)dX{u) 


( 6 . 6 ) 


so that 



f{v - w) g{-w) dA(w;) dy^^^{v) 


{f *g){v)dy^^^{v) = ylii{f*g), 


( 2 ) 


Z^leU/*/) = = E[//d^ J/d^ 


= E 


ffd<P ffd<P 


= E[|^(/)|2] >0. 


), 


Remark 6.1 (see also [H Rem. 13]). One can alternatively define 

X A) = E[¥{A)(P{A')] for bounded A, A G 
and then obtain /ired^'^ from this as above, via 

/ f{x,y)dy^‘^'^^^\x,y) = [ [ f {u + v, u) dX{u) 

f{u,u-v) dA(M). 
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Then, we have and = /^red- Since 


,( 2 ) 




f{y,x)dfx^'^\x,y) = / f{x,y)dy(^'>{x,y) = 


f{x,y) dy^'^’^^^\x,y), 


we see that the alternative choice of factoring out the shift invariance in Eq. fl6.4p . namely 
integrating f{u,u + v) on the right-hand side of this equation, leads to where 


f(x,y)Aiif^\x,y) = / / /(ti, u + ii) dA(ti) 


(6.7) 


We choose the dehnitions as in Eqs. (Ih.dh and 


because these fit well to the formulation 

(2) ,,(2.alt) 


of the limit in Eq. fIb.Sp below. Note that, in the real-valued case, and agree. 


The ‘complex-valued’ analogue of [U Thm. 5] now reads as follows. 


Theorem 6.2. Assume that conditions fl 6 .ip and fl 6 . 2 p are satisfied, and let (Pn '■= 
denote the restriction of <P to the open ball of radius n around 0. Then, the natural 
autocorrelation of T>, which is defined with an averaging sequence of nested, centred balls, 
almost surely exists and satisfies 


( 6 . 8 ) 


(<p) 

7 ^ W= lim 

n^oo XyBfi) 


where the limit refers to the vague topology, 
random. 


lim = 

n^oo 

In particular. 


Tved'! 

the autocorrelation is non- 


Proof. The proof is a suitable ‘complex-valued interpretation’ of the proof of [U Thm. 5]. 
Fix a continuous function /: M'’* —> C with compact support. We have to check that 

(6.9) 7^ d2,(/) (a.s.). 

Since both sides are locally hnite (complex-valued) measures, it actually suffices to check 
Eq. fl6.9p for real-valued /. For x G M'’*, dehne 


F{x) := [ f{x-y)d^{y) = [ f{x + y)d^{y). 

Jra J^d 

Clearly, F inherits stationarity and ergodicity from wherefore F is a (complex-valued) 
ergodic random function on in the sense that shift-invariant events for F have ‘trivial’ 
probabilities (0 or 1 ), and we obtain 


E 



F{x)\d\d>\{x) 


< oo 


for any bounded and measurable A C M'’*. 
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Define a (complex-valued) additive covariant spatial process X{A) in the sense of [66] , 
indexed by a bounded and measurable A C via 

X{A) := f F{x)d<P{x). 

Ja 


Covariant in this context means that X behaves ‘naturally’ under translations: When R'’* 
acts on X via {TuX){A) := f^F(x) d(T„^)(a;), for u G R'’*, then (TuX)(A + u) = X{A). 

Decomposing X into its real and imaginary parts (by decomposing F and ^ and suitably 
grouping terms) we can apply [661 Cor. 4.9] to obtain a.s. 


lim 

n^oo 


= lim 


A(5, 


F{x) d^{x) 


^ Bn 


lim 

n^oo 


X{B^) 

KBn) 


= E 


X^ 

A(5i) . 



f{x - y) d^{y) d^(x) 


1 

Wi) 

1 

Wi) 


Isi {x)f{x-y) d^^‘^\x,y) 



{x)f{z) d/iJ^W dA(x) 



The difference between ^ and <Pn * is a (random) ‘boundary term’ that almost 
surely vanishes in the limit as n —)■ cxo. To prove this formally, decompose ^ -|- i^^., 

— i^Q and then argue as in the proof of [H Thm. 3] for each of the four terms 
appearing in ■ D 


Remark 6.3. Theorem 16.21 allows to reformulate Theorem 4 and Corollary 1 from |1] for 
complex-valued clusters as follows. If ^ is a stationary ergodic point process, i.e. ^ is a 
random sum of Dirac measures, with distribution P satisfying Eq. fl6.2p . and if we replace 
each point independently by a random complex-valued measure with distribution Q, then 
the formulas describing the autocorrelation and the diffraction of the resulting cluster 
process given in [H Thm. 4 and Cor. 1] continue to hold. 


Let us also mention that, by specialising ^ to a renewal process. Theorem 16.21 allows to 
recover Eq. fl4.ip and, in particular. Theorem 14.11 from this more general perspective; see 
[1] for further details, and how this can be used to formulate the renewal process also for 
more general ‘dropping’ distributions. 


6.1. A ‘Palm-type Distribution’ for Complex-valued Random Measures. In the 

case of a positive random measure Eq. (16.Sp can be interpreted via the Palm distribution 
Pq of the law of which is a probability measure on locally finite measures (intuitively, 
the law of ^ viewed relative to a typical point of its support) via 

(6.10) = pJpo 
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where p > 0 is the intensity and ^Po the hrst moment measure of Pq! compare [U Sec. 5.2], 
This interpretation breaks down in general in the signed or complex-valued case because 

f2) 

/ijigj will not be a positive measure. One way to extend this line of thought is to re-interpret 
the Palm distributions in a way suited for complex-valued random measures as follows. 

Recalling the structure of the polar decomposition, the random measure ^ can equiva¬ 
lently be described via (|^|, ^ph), where |^| is the total variation measure and the mapping 
(Pphi —)■ [0,27r) the ‘phase function’: 

(6.11) f f{x)d<P{x) = f f (x) d\<P\{x). 

Note that $ i—)■ (|^|,^ph) is measurable, so (|^|,^ph) is in fact a random variable. Define 
a positive cr-hnite measure C on M'’* x (this is the equivalent of the so-called Campbell 
measure for the complex-valued context and agrees with the usual Campbell measure if 
is a positive random measure) via 


g{x,Lp)dC{x,Lp) := E 




g{x,e d|^|(a:) , 


whenever the right-hand side is defined (which will for instance always be the case when g 
is measurable and non-negative). By the shift invariance of and hence that of |^|, the 
projection of C to M'’* is p times Lebesgue measure (with p G [0, oo) being the intensity of 
|(P|), hence there is a family of probability measures Px on Al, with Pj, G V{M.) for all 
X G M'’*, so that we can disintegrate (compare [ISl Thm. 15.3.3]) 


( 6 . 12 ) 


gdC = 


g{x,ip) dPx{(p) pdX{x). 


J JRrf Jp{M) 

Definition 6.4. We call the elements of the family {Pj, | a; G M'’*} the Palm distributions 
in the complex-valued case. 

Let, for A C bounded and measurable. 


IpM) ■■= [ AA)dPx{^) 

Jm 

be the expectation (or hrst moment) measure of Px- By shift invariance, we have Pj, = TxPq, 
X G M'’*, and hence Ip^ = Txlp^. The connection between the (reduced) second moment 
measure and the Palm distribution carries over to the complex-valued case as follows. 


Proposition 6.5. For the extended definition of the Palm distribution, one has 

(2) _ r 
Fred — P^PqI 


SO Eq. fl6.10l] also holds in this case. 
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Sketch of Proof. Consider g{x, ip) = lA'{x)ip{A) with 1 denoting the characteristic fnnction 
and with A, A' C bonnded and measnrable. Then, 


( 6 . 13 ) 


(yfdC = E 


lA'ix)e ^{A) d\<P\{x) 


= E 


^{A) / lA/(x)e-‘‘^p>'(^)d|^|(x)j = E[<P{A)^{A')] 


by dehnition, whereas the disintegration formnla yields 


(6.14) 


gdC= / 1a'{x)p>{A) dPx{p)) p d\{x) =P Ip^{A)dX{x) 

jRd Jm JA' 


= P Ipa{A - x) dX{x) = p 


'A' 



lA'{x)lA-x{y) d/po(2/) dA(a:) 


= P 


= P 




lA'{x)lA{y + x) dA(a;) dlp^{y) 

lAi-x)lA{y - x)dX{x)dIp^{y) = p Ipq{1a * 1a') ■ 


Comparing Eqs. fl6.13p - fl6.14p with Eq. fl6.6p yields the claim. 


□ 


If ^ is ergodic, the viewpoint that Pq describes the conhgnration relative to a point in 
the snpport drawn according to ^ is corroborated by 

d\<P\{x) f g{p,)dPo{p) (a.s.) 

A[.Dn) Jb„ Jm 

for any bonnded measnrable g: A4 —h M. 

The viewpoint of (possibly complex-valued) ergodic random measures for diffraction is 
a useful one since it provides a connection to the large literature on random measures 
and on stochastic geometry; see [30], HO] [5T1 (52] [62l H] and references therein, as well as 
ua for a recent generalisation that can also be considered from the diffraction point of 
view. However, our approach also shows a limitation that one encounters when trying to 
infer properties of a random conhgnration of scatterers from its kinematic diffraction: As 
is evident from Eq. fl6.8l) in Theorem 16.21 the only ‘datum’ from a random ^ visible in 
its autocorrelation, and hence also in the corresponding diffraction, is the second moment 
measure. It is well known that second moments are generally insufficient to determine 
the distribution of ^ unless further structural properties are known. This inverse problem 
is known as the homometry problem in crystallography and the inference problem in the 
theory of stochastic processes. 
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7. Outlook 

Our exposition provides a snapshot of the present knowledge about systems with con¬ 
tinuous diffraction components; see [121 Chs. 10 and 11] as well as [H [15] for additional 
examples, and [sniiisi for connections with the dynamical spectrum. Nevertheless, as is 
apparent from a comparison with the pure point diffraction case [2311131 ISSl EH EH], the 
status of general results is lagging behind. Even for many important examples, some of 
the most obvious questions are still open from a mathematical point of view. In particular, 
this is so for random tiling ensembles in dimensions d > 2, or for equilibrium systems just 
beyond the complexity of the (planar) Ising model. 

Apart from the systems considered here, an interesting class is provided by random 
substitution and inflation systems, as introduced in [JT]. The randomness present here is 
compatible with the long-range order of Meyer sets with entropy [211 IMl ES] , which means 
that one obtains interesting mixtures of pure point and absolutely continuous diffrac¬ 
tion measures. Though this direction has not attracted much attention so far, it is both 
tractable and practically relevant. 

From a more general perspective, one lacks some kind of analogue to the key theorems in 
pure point diffraction (such as the Poisson summation formula or the Halmos-von Neumann 
theorem). While there is at least the theory of Riesz products [231 ES] for self-similar 
systems with singular spectra, a general approach to stochastic systems is only at its 
beginning. Methods from point process theory [30], such as the Palm measure and its 
connection to the autocorrelation (via its intensity measure), look promising, but have 
not produced many concrete results so far. The latter, however, are needed to make 
some progress with the complicated inverse problem for such systems. Though there is 
substantial knowledge from the inference approach [19], it is not clear at present how this 
can be used, and how reasonable restrictions could be included. 
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